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Abstract 

Several astrophysically relevant solid oxides and silicates have extremely small opacities 
in the visual and near-infrared in their pure forms. Datasets for the opacities and for 
the imaginary part k of their complex indices of refraction are hardly available in these 
wavelength ranges. We aimed at determining k for spinel, rutile, anatase, and olivine, 
especially in the near-infrared region. Our measurements were made with impurity- 
containing, natural, and synthetic Stardust analogs. Two experimental methods were 
used: preparing small sections of natural minerals and synthesizing melt droplets under 
• the electric arc furnace. In both cases, the aborption properties of the samples were 

measured by transmission spectroscopy. For spinel (MgAl204), anatase, rutile (both 



m 

^- 

Q ■ Ti02), and olivine ((Mg,Fe)2Si04), the optical constants have been extended to the visual 

^ I and near-infrared. We highlight that the individual values of fc(A) and the absorption 

^ i cross section QabsW depend strongly on the content in transition metals like iron. Based 

on our measurements, we infer that k values below 10~^ are very rare in natural minerals 
including Stardust grains, if they occur at all. Data for k and Qabs{^) are important for 
various physical properties of Stardust grains such as temperature and radiation pressure. 
With increasing Qa6s(A) due to impurities, the equilibrium temperature of small grains 
' in circumstellar shells increases as well. We discuss why and to what extent this is the 

\^ ■ case. 



1 Introduction 

Most studies of the absorption properties of Stardust analogs are focused on those wave- 
length ranges where the respective minerals or glasses have strong resonance features. 
This is usually the case in the ultraviolet range, on the one hand, and in the mid-infrared 
(MIR) range, on the other. Between these two ranges - i.e. between the range of strong 
electronic resonances and the range of strong lattice vibration resonances - there is fre- 
quently a lack of data for the complex index of refraction m(A) = n(A) -I- i fc(A). This is 
especially true for those insulator materials that are highly transparent in the visible and 
near infrared (NIR), such as magnesium and aluminum silicates and oxides, since their 
/c(A)-values tend towards values that are too small to be measured with respect to several 
standard methods of spectroscopy. For example, the fc(A)-values of pure MgAl204-spinel 
would decrease from 10"'^ at about 7fj,m to 2.2 x 10"^'^ at about 2/im according to Tropf 
and Thomas (1991). For A < 2/im, no /c-data are available at all from this source. 

The importance of A:(A)-data characterizing Stardust in the visual and NIR region 
becomes clear by the following argument. Dust in circumstellar shells mostly absorbs ul- 
traviolet, visual, and NIR radiation from the stellar radiation field, thermalizes it via its 
internal degrees of freedom and re-emits it at mid-, as well at far-infrared wavelengths. 
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Consider an asymptotic giant branch (AGB) star with an effective photospheric tem- 
perature of 2500-4000 K: such an object will emit most of its radiation energy close to 
1 /Ltm. Having no information on k{\) for A « 1 /xm means nothing else than being unable 
to calculate the energy budget of a dust-enriched stellar atmosphere. Therefore, efforts 
should be made to close the gaps in A: (A) - which are usually much wider than eventual 
gaps in n(A) - where they show up in optical databases for Stardust components. 

For a very limited number of potential cosmic dust minerals - namely for spinel, rutile, 
anatase, and olivine - this task will be approached in the present paper by measuring the 
absorption properties depending on their content of transition metal ions as 'impurities' 
(colouring and absorption-enhancing). 

Why do standard optical databases - especially those used in astronomical modeling 
- often lack absorption data for the transparency regions? Reflectance spectroscopy on 
bulk surfaces, which is very helpful for deriving MIR optical constants, hardly deliver 
any reliable constraints on A; (A) in the transparency regime. Here, absorption coeffi- 
cients based on reflection measurements are usually nothing more than the results of 
extrapolations. A proper determination of the usually small absorption coefficients in 
the near-infrared range require transmission measurements at high column densities. 
Since scattering losses may compete with the eff'ects of absorption in the visual and NIR 
region, homogeneous (bulk) samples are needed, i.e. sections of millimeter or tenths of 
millimeter thickness. These are not easy to obtain whether synthetically or from natural 
sources. 

This paper intends to provide a guideline of how to improve A; (A) data for the task 
of modeling dust absorption in the visual and NIR wavelength regions. The paper is 
structured in the following way. Section 2 aims at characterizing some of the absorption 

processes relevant to these regions. In Section 3, we describe our measurements on 
individual Stardust analogs, namely spinel, rutile, anatase, and olivine. We present the 
spectra and optical constants derived for these materials, focusing ou their transparency 
regimes. In the final Section, we disciiss a specific astrophysical implication of our work, 
the influence of the derived A;(A)-values on dust temperatures in a radiative equilibrium. 



2 The concept of the 'transparency regime' 
2.1 Definition and illustration 

Figure 1 shows k{X) for glassy Si02 according to Philipp (1985), for Mgo.5Feo.5O accord- 
ing to Henning et al. (1995), and for MgO according to Roessler and Huffman (1991), 
from the ultraviolet to the far infrared (FIR). Between the two ranges of strong absorp- 
tion, where k{\) reaches and exceeds unity, a large region in the visual and near-infrared 
(NIR) can be seen, where fc(A) drops to very low values both for Si02 and MgO, while 
this is not the case for the iron-rich material Mgo.5Feo.5O. Recalling that the absorption 
coefficient a(A) of a solid, which enters into Lambert-Beer's law, is 

47rfc(A) 

a(A) = — - — , (1) 

we can quantitavely define, following Mitra (1985), the 'transparency regime' of a crys- 
talline or amorphous solid as the region where 

a{X)<lcm-\ (2) 

At a thickness of 1 mm, a window in this regime would transmit more than about 
90% of the light entering at normal incidence; at 1 cm thickness, it would still transmit at 
least 1/e or 37% of the radiation (neglecting reflection losses). In terms of A;(A), eq. (2) 
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Figure 1: Absorption indices k{X) for glassy Si02, crystalline Mgo.5Feo.5O, and crystalline 
MgO from the ultraviolet to the far infrared. The grey area illustrates the 'transparency 
regime' as defined by Mitra (1985) and by equation (2). The sources of the data sets are: 
[1] Philipp (1985), [2] Henning et al. (1995), and [3] Roessler & Huffman (1991). 



defines a wavelength-dependent limit as shown in Fig. 1 (border between the gray and 
white areas). At A = 1 /zm, for example, we may call a solid 'transparent' if fc <8 x 10~^. 
This is definitely the case for pure MgO and Si02. It is also true for pure crystalline 
MgAl204, a-Al2 03, Fe-free Mg-silicates and indeed most oxygen- rich Stardust analogs 
without impurities. We note that for the visible/near-infrared, our transparency criterion 
is consistent with the statement by Bohrcn & Huffman (1983) who point out that already 
for fc-values amounting to kyis « 10^^ bulk materials appear black, even in small pieces. 

While the transparency regime is limited towards shorter wavelengths in the UV/ visible 
by the edge of electronic interband transitions and possibly excitonic transitions (we give 
the energies and wavelengths at which strong absorption sets in Table 1), in the infrared 
it is limited by the energetically highest lattice vibration mode. As a rule of thumb, if a 
solid has a transparency region, it sets in at wavelengths shorter than 

A* = 1/2Alo = 1/(2wlo), (3) 

where ojlo denotes the highest longitudinal optical frequency (in wavenumbers) and Xlo 
denotes the corresponding wavelength (cf. Barker 1975; however, Barker uses a different 
criterion for the transparency and therefore derives other values for A*). Physically, 
this upper frequency limit corresponds to the limit up to which two-phonon absorption 
processes can be excited. Photons of shorter wavelengths (i.e. in the transparency region) 
can be absorbed by higher overtones, i.e. the simultaneous excitation of three or more 
phonons, which is less likely and does usually not lead to k values above the transparency 
limit. As can be seen from Table 1, transparent materials like MgO, amorphous Si02, 
MgAl204, and Ti02 have two properties in common: their highest energy LO mode 
is located at 700 to 1400 cm~^, and accordingly, their transparency regime sets in at 
wavelengths below 7.1 to 3.6 /im. This is also the wavelength region where many sets of 
optical constants derived from infrared measurements come to an end because there is a 
lack of physically meaningful k values (as mentioned above) . 

For nonconducting materials, there is an FIR transparency regime as well. At wave- 
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lengths greater than a few hundred micrometers, k again drops to values lower than 
10~^ to 10~^, which are low enough to fulfill (2). Measurements in this wavelength 
range, therefore, bear some similarities to near-infrared measurements, although wave- 
lengths and absorption mechanisms are different. Data in this wavelength range are 
additionally limited by the availability of FIR-spectrometers. 

2.2 Causes of transparency and of absorption in the NIR region 

It can furthermore be seen from Fig. 1 that there are materials without any transparency 
regime within the wavelength range shown. Mgo.5Feo.5O a solid solution of MgO and 
FeO - is an example of such a material. FeO is a semiconductor, in contrast to MgO and 
Si02. Actually, it is a nonstoichiometric compound, its sum formula should better be 
written as Fci_2;0. Therefore, it contains both Fc^+ and Fe^+ ions, the content of Fc'^+ is 
0.09 according to Henning et al. (1995). These two properties lead to various absorption 
mechanisms for wavelengths in the visible and near infrared, as well as the FIR, which 
arc not present or occur at different wavelengths for the previously discussed iron-free 
oxides, but are efficient in all the members of the series Mga;Fei_a;0 (x<l). Among these 
mechanisms are 

1. the fundamental interband absorption edge occurring in the visible wavelength 
range (band gap 2.3 eV for FeO, according to Park et al. 1999); 

2. free electron excitations causing the increase of fc(A) towards the FIR. For FeO, 
this absorption mechanism should be weak, but a part of the material, which is 
metastable at room temperature, may be decomposed into Fe and Fe304, which 
both have high free-electron densities (at room temperature for Fe304); 

3. charge transfer transitions between the orbitals of Fe^+ and those of other metal 
or ligand ions, that are an important source of strong absorption in the visible 
extending into the near infrared; 

4. the splitting of the d-electron energy levels in transition metal ions, depending on 
their coordination by oxygen ions. In case of lower symmetry, transitions between 
these d-clcctron states can lead to moderately strong absorption bands (so-called 
crystal field bands. Burns 1993). 

The last three mechanisms can obviously enhance the absorption in otherwise trans- 
parent wavelength ranges above the transparency limit. All three of them and even more 
absorption mechanisms, such as vibrational excitation, can be introduced by impurities. 

If a (stardust) mineral is transparent at near-infrared and visual wavelengths in its 
pure, ideal, stoichiometric form, it may in reality contain 'impurities' that strongly in- 
crease fc(A). This phenomenon is well known in geochemistry and mineralogy, but rarely 
accounted for in astromincralogical applications. The most efficient 'enhancers' of fc(A), 
consequently of mass absorption coefficients (opacities), are transition metals: Sc, Ti, 
V, Cr, Mn, Fe, Co, and Ni. Among them, Fe, Ni, Cr, and Mn are probably the most 
relevant for stardust minerals considering the mean cosmic elemental abundances (Burns 
1993). 

Any significant inclusion of transition metal ions into an originally transparent mineral 
will turn its appearance in the visual spectral range from colorless to colored or even 
black. At present, it is hardly known which amounts of transition metal ions are actually 
present in stardust minerals as 'impurities', but it is unrealistic to expect impurity-free 
minerals in circumstellar and interstellar environments. 

Even though several visually transparent oxides such as spinel and AI2O3 have been 
identified in presolar grains, it has not yet been possible to determine the transition 
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Table 1: Electronic band edges including strong excitons in terms of energy Eg, wave- 
length Xg, longitudial optical lattice frequencies cjlo, and onset wavelengths A* for the 
transparency regions for different materials (see text for more details). 



Material 




A. 


UJLO 


A* 


Ref. 




[eV] 


[/im] 


[cm^-'-] 


[^m] 




Am. Si02 


8.3 


0.15 


1245 


4.01 


[1] 


MgO 


7.5 


0.17 


728 


6.87 


[2] 


MgAl204 


7.75 


0.16 


877.2 


5.7 


[3] 


Mg2Si04 


7.5 


0.17 


1078 


4.6 


[4],[5] 


Ti02 rutile 


2.9 


0.43 


831.3 


6 


[6] 


Ti02 anatase 


3.2 


0.39 


872.8 


5.7 


[7],[8] 


FcO 


2.3 


0.54 


526 


9.5 


[9], [10] 



[1] Philipp (1985); [2] Roessler & Huffman (1991); [3] Tropf and Thomas (1991); [4] 
Shankland (1968); [5] Sogawa et al. (2006); [6] Ribarsky (1985); [7] Tang et al. (1993); 
[8] Posch et al. (2003); [9] Park et al. (1999); [10] Henning et al. (1995). 



Table 2: Mean cosmic elemental abundances of selected transition metals according to 
Palme and Beer (1993). 



Element 


abundance log N 




(N(H)=10i2) 


Sc 


3.1 


Ti 


4.96 


V 


4.00 


Cr 


5.68 


Mn 


5.45 


Fe 


7.49 


Co 


4.91 


Ni 


6.24 
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metal content of these minerals; in other words, their 'color' has not been quantitatively 
measured. Spectroscopic methods determining the transition metal content of cosmic 
dust by 'remote sensing' are restricted, up to now, to reflectance spectroscopy of lunar 
and planetary surfaces (Burns 1993). 

Therefore, in the present paper, we chose the following approach when synthetically 
including transition metals into oxides: We based our measurements of k{X) partly on 
natural terrestrial samples, which always contain impurities to some extent, and we were 
partly able to widen the range of possible inclusions of transition metals by producing 
synthetic samples. For one of our samples, spinel, it was thus possible to create a series of 
samples and corresponding sets of optical constants such as exploring the possible range 
of absorption that may occur in samples of cosmic elemental composition. 

3 Individual measured Stardust analogs 

3.1 Derivation of k from the transmission spectra 

The method we applied to derive A;(A) was transmission spectroscopy of small sections 

of our materials. For producing these sections we embedded pieces of our materials with 
sufficient size (diameter > 0.5 mm) into resin and cut and abraded the resin slabs to 
thicknesses between 80 and 1100/um. After that, the samples were just polished to have 
an even and smooth surface. 

Some of the examined materials are crystals of non-cubic symmetry and therefore 
optically anisotropic; i.e. their optical constants depend on their orientation relative 
to the polarization direction of the light. For these materials, it has been necessary to 
perform measurements with polarized radiation. An overview of the sample properties 
such as chemical composition determined by energy-dispersive X-ray analysis (EDX), 
crystal symmetry, and slab thickness is given in Table 3 (the oxygen content is not a direct 
result of the EDX measurements, but has been inserted according to stoichiometry). 



Table 3: Compositions of the samples and thicknesses of the platelets examined in this 
paper. 



sample 


chemical composition 


thickness 


crystal symmetry 


Natural Spinel^ 


Mgi.o2Ali.93Cro.o087Feo.oi204 


250/im 


cubic 


Synth. Spinel 1 


Mg1.00Al1.98Cr0.02O4 


80/im 


cubic 


Synth. Spinel 2 


Mgo.93Al1.99Cro.03O4 


270/xm 


cubic 


Synth. Spinel 3 


Mgo.98Al1.89Cro.12O4 


120/im 


cubic 


Synth. Spinel 4 


Mgo.98Al1.78Cro.23O4 


BO^m 


cubic 


Natural Rutile^ 


Tio.984Vo.o08Feo.o08 02 


155/;im 


tetragonal 


Natural Anatase 


Tio.992Vo.o0802 


255 iim 


tetragonal 


Olivine San Carlos 


Mgl.96Feo.16sio.8904 


1063/1115Mm 


orthorhombic 


Olivine Sri Lanka 


Mgl.56Feo.4sio.9104 


300/xm 


orthorhombic 



f The natural spinel and rutile samples have small (< 1%) OH contents that are not 
shown in the formula. 



For the spectroscopic measurements, the following spectrometers have been used: a 
Perkin Elmer Lambda 19 (for the wavelength range from 0.275 to 2 fJ,m), a Brukcr FTIR 
113v (for the wavelength range from 1.6 to 25 /^m), and a single-beam optical absorption 
spectrometer constructed on the basis of a SpectraPro-275 triple-grating monochromator. 
The actual range of the transmission measurements is limited by the wavelengths at which 
the samples become opaque. For a typical sample thickness of 300 /xm, this happens for 
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a k value of about 10"'^ on the short-wavelength side and at a few 10~^ on the long- 
wavelength side. The transmittance values (transmitted intensity normalized by incident 
intensity) were converted in k{X) using the relation 

where d is the thickness of the small section and the reflectance R is given as 

_ (n(A) - 1)^ + fc(A)^ 

KA) + l)2 + fc(A)2- ^ ' 

For very small values of k(X) this equation reduces to 

^-Ra)TTF- 

While using ri(A) values from Palik (1985) we were able to derive values for fc(A) for 
various samples. 

The following subsections contain a presentation and discussion of our results for the 
individual materials. The derived k values for the NIR will be made publicly available on 



the Jena Database of Optical Constants, http: / /www. astro. uni-jena.de/Laboratory/Database/databases. html 



3.2 Spinels 

3.2.1 Material properties, sample preparation and analysis 

Spinel is an abundant component of presolar grains from meteorites (e.g. Hoppe 2004) and 
has been proposed as the carrier of mid-infrared emission bands observed in the spectra of 
circumstellar dust shells (e.g. Posch et al. 1999, Fabian et al. 2001). Pure stoichiometric 
spinel (MgAl204) has an elementary cell that consists of 32 anions, 16 Al, and 8Mg 
cations. In the ideal crystal lattice, the trivalent Al cations are octahedrally coordinated 
by oxygen ions, while the bivalent Mg cations are tetrahedrally coordinated. The cubic 
lattice symmetry of this mineral make it optically isotropic. 

While many previous publications refer to impurity-free spinels, the measurements 
presented in this paper refer to a series of spinels that contain impurities (mainly chromium 
and iron). The Cr ions partly substitute Al, so they occupy a part of the octahedral sites. 
The Fe ions partly substitute Al, partly Mg. Because the lattice structure of the Cr- 
containing spinels is unambiguous (since Cr is not polyvalent in contrast to Fe), we have 
chosen synthetic Cr-doped spinels as a representative of the effect of impurities on k, even 
though in cosmic environments, it is much more likely to have Fe impurities in spinels 
due to the much higher stellar abundance of Fe compared to Cr as given in Table 2. 

Additionally, a dark-red natural spinel crystal from Burma was examined. Apart 
from a very small amount of Si, it contains traces of both Fe and Cr. The exact chemical 
composition of the sample is given in Table 3. 

The synthesis of the Cr-doped spinels was performed in the following way (see also 
Richter et al. 2005). We used a tungsten electric-arc furnace with a water-cooled cop- 
per baseplate. Powders of Mg oxide, Al oxide, and Cr(III) oxide have been mixed in 
corresponding proportions so as to obtain the basic material for the melting procedure. 
After being homogenized, the mixture was pressed into several pellets of 150-300 mg 
mass each and a diameter of 1.3 cm each. This densification of the powder was neces- 
sary to avoid electrostatic interactions between the tungsten pike in the furnace and the 
powder particles when starting the melting process. Since Mg oxide in particular is a 
hardly fusable material, the pellets were broken again into several pieces (4-7 for each 
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Figure 2: Section of the spinel crystal from Burma used for deriving k{X) from transmis- 
sion spectroscopy. From bottom to top this crystal has a size of around 3 mm. 

pellet) to facilitate the melting process. After putting several pieces into the furnace, the 
melting process started under a 1 bar argon atmosphere. The samples were exposed to 
the electric arc until they were fully melted. The formation temperature should be close 
to the melting point of spinel (about 2400 K), and it probably increased with the Cr 
content (the melting point of Cr203 is about 2700 K). After cooling to subliquidous tem- 
perature within a few seconds and subsequent cooling to room temperature, we obtained 
melt-droplets that were 3-5 mm in size and had colors ranging from light red in the most 
Cr-poor spinel to dark red (almost black) in the most Cr-rich sample. The theoretical 
compositions of the samples are given by MgAl(2_2;)Cr2;04, with x amounting to 0.02, 
0.03, 0.12, and 0.23. 

The exact chemical compositions were derived from EDX analyses, where several 
(typically 12) spots for each of the Cr-spinel droplets were scanned and investigated for 
inhomogeneities and variations in the Mg/Al/Cr ratios. Tab. 3 gives the resulting mean 
stoichiometries - also for the Ti02 and olivine samples discussed in Sections 3.3 and 3.4. 

3.2.2 Measured transmission spectra and derived k values 

We now present and interpret the results of the transmission spectroscopy performed 
on the above-mentioned samples at UV-, visual, near-infrared, and mid-infrared wave- 
lengths. 

The measured transmittances were converted into k{X) values by applying (4) and 
Fig. 3 shows the derived results. The dashed line represents k{X) for pure stoichiometric 
spinel according to Tropf and Thomas (1991). The mid infrared maxima of k(X) - 
between 10 and 40/im - correspond to three vibrational modes of the lattice structure: 
The origin of the 31-32/im band is an Mg-0 vibration, while the remaining two maxima 
originate from Al-0 vibrations. For further details see Fabian et al. (2001). 

It can be clearly seen from Fig. 3 that for the impurity-free spinel - and only for this 
one - a very steep decrease of k{X) with decreasing wavelength occurs at A < 7/im. 

The solid line denotes k(X) for our natural sample from Burma which contains about 
1% iron and approximately as much chromium. This A:(A)-curve begins to deviate from 
the corresponding line for pure spinel at 6-7 ^m and never drops below the limit of 
the transparency regime, in accordance with the opaque character of the 0.25mm thick 
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Figure 3: Absorption indices A; (A) of different spinels derived from transmission spec- 
troscopy. For the interpretation of the different bands seen in the spectra, see text. 

sample. Within the range of rough agreement between the two curves (> 7/im), there is 
some structure in our new data, especially a peak at ~7 fxm, which is likely an overtone 
of the 13-14 /nm lattice vibration band. 

Between 4 and 1.5 /xm, a broad maximum occurs with its origin in the Fe^+ and 0H~ 
content of the sample. Even though the OH content is very small (definitely below 1%), 
it is sufficient for the 0-H stretching vibration band to appear in the spectrum as a 
narrow band at 2.8 /xm. Then fc(A) approaches values even close to 10~^ at the center of 
this band. The same holds true for the maximum of the Fe^+ band, which is located at 
2.0 fim. A detailed interpretation of the Fe^^ and 0H~ bands in natural and synthetic 
spinels can be found in a recent paper by Lenaz et al. (2008). A comparison between our 
data and those published by Lenaz et al. (2008) is shown in Fig. 4. 

The region below 1 /um is characterized by further electronic absorption mechanisms. 
In this region, the fc(A) spectra strongly depend on the Cr content of spinel. This can 
be nicely seen by comparing the solid line with the dotted lines, which represent the 
Cr-doped synthetic spinels. 

The Al to Cr substitution in the samples results in a generally increased visual and 
UV absorption. Additionally, two distinct absorption bands occur around 0.550 /im and 
0.387 /xm. The positions and shapes of these bands accordance well with the band posi- 
tions in spectra of low-Cr Mg-Al- spinels - both natural and synthetic ones - published 
by Wood et al. (1968), Taran et al. (1994), and Ikeda et al. (1997). These bands originate 
from crystal field electronic transitions, namely '^T2g '^■^2g for the 0.550 /Ltm maximum 
and *Tig *A 2g for the 0.387 /um maximum. Of course, the intensity of the Cr elec- 
tronic absorption bands is higher for the synthetic Cr-doped samples, and fc-values up to 
2 X 10""^ are reached at 0.550 /xm for 12% Cr content. 

3.3 Rutile and anatase 

Together with Ca and Al oxides, Ti oxides are considered as candidates for the very first 
condensates in oxygen-rich circumstellar shells, notwithstanding that titanium is rather 
rare according to its mean cosmic elemental abundance (2.5 orders of magnitude less 
abundant than Si and Mg). With respect to dust formation in circumstellar shells, it is 
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Figure 4: Comparison between the absorption indices derived for spinels with different 
transtion metal content. The solid line refers to our sample - with an Fe content of 1.2% 
and a Cr content of 0.87%. The dotted line refers to a very Fe-rich sample (Fe-content 
of 11.1%), the dashed line denotes an Fe-poor sample (both according to Lenaz et al 
(2008)). The Fe content of the samples clearly correlates with the respective fc- values at 
2/xm. 



noteworthy that molecular TiO is extremely prominent in the atmospheres of M- and 
S-stars. Furthermore, titanium is heavily depleted in the interstellar gas, indicating that 
it is indeed consumed by previous dust formation processes. 

Jeong et al. (2003) predict that Ti02 is the most promising candidate for the first 
(seed-)condensate in oxygen-rich circumstellar shells. There is also evidence from me- 
teoritics that Ti compounds form in circumstellar shells and a presolar meteoritic Ti02 
grain has been tentatively identified (Nittler 2003). 

We examined the two most common terrestrial species of Ti02 with respect to their 
near infrared absorption properties, namely anatase and rutilc. Both materials have a 
tetragonal crystal symmetry, hence only one optical axis. Consequently, the measure- 
ments for the polarization directions parallel and perpendicular to this optical axis could 
be performed on a single sample cut along this axis by rotating the polarizer by 90 
degrees. 

k{X)NiR is not unmeasurably small in the case of our samples due to their Fe and 
V contents. The Fe- and V- contents of our Ti02 samples have been determined by 
energy-dispersive X-ray measurements and are summarized in Table 3 



3.3.1 Rutile 

Rutile is the most abundant modification of titanium dioxide on the Earth. It has 
tetragonal symmetry like anatase, but is characterized by a higher density, a greater 
Ti-Ti distance, and a shorter Ti-0 distance than the latter. As for rutile's ultraviolet 
and mid-infrared properties, Ribarsky (1985) compiled a rather comprehensive table of 
optical constants. However, for wavelengths between 0.4 /xm and 11.2 /xm, data for A;(A) 
are largely missing. 
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Figure 5: Section of the rutile crystal used for deriving k{X) from transmission spec- 
troscopy. The crystal has a width of about 1.5 mm. 

We were able to measure k for a natural rutile sample from 0.5 to 8 fiui, based 
on a section with a thickness of d=155/im (see Figs. 5 and 6). The most significant 
feature in the NIR fc-spectrum of rutile is a sharp peak at 3.03 fj,m. This peak, which 
becomes especially prominent when the electric field vector of the polarized radiation is 
perpendicular to the c-axis of the crystal, is due to an OH stretching mode (see Maldener 
et al. 2001). It is explained by the presence of H2O in natural rutile. The 3.03 /im band 
of our sample may be explained with an H2O content of the order of magnitude of 300 
weight ppm, based on the spectra and analytical results by Maldener et al. (2001). 

3.3.2 Anatase 

Anatase is a modification of tetragonal Ti oxide with a TiOe octahedron structure con- 
sisting of one Ti'*+ and six O^^ ions. It can be transformed to rutile by heating above 
1200 K. Its mid-infrared optical constants have been published by Posch et al. (2003). 

For anatase we used a section of d=255 fim thickness to perform transmittance mea- 
surements between 0.7 and 8/im. The transmittance spectra obtained for the individual 
orientations of the sample relative to the polarized radiation did not differ within the 
measurement's accuracy, therefore only one curve is plotted in Fig. 7. 

The most conspicuous feature in the fc(A) spectrum of anatase is the broad double- 
peaked band 2-3 /im. Maximum k values of more than 0.001 are reached in this range. 
The second of the two peaks, located at 3.0 /im, corresponds to an O-H stretching mode 
at 3389.7cm— 1. It is an indirect consequence of the V content of our sample, since V- 
containing Ti02 has a stronger hydroxyl content (Zhou et al. 2010). The primary peak 
in the absorption, located at 2.0 /im, may be directly caused by the V content of the 
sample. 

3.4 Olivine 

Olivine is, together with pyroxene, the most abundant crystalline silicate in space. Its 
prominent bands are observed in many spectra of accretion disks around young stars 
and in the outflows of AGB stars (Henning 2003a, 2003b). Olivine also occurs in our 
solar system in the form of cometary dust (Hanner 2003), as interplanetary dust particles 
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Figure 6: Absorption index A:(A) for rutile (Ti02) derived from transmission measure- 
ments on a small section of 155 /im diameter. The fc(A) values in the UV and MIR are 
from Palik (1985-98). The sharp peak in fc(A) comes from a small crystal water content 
of the natural sample. 



present paper 

k parallel, Jena database 

k perpendicular, Jena database 

grey; transparency regime 





X[Mm] 



Figure 7: Absorption index k{X) of the Ti02 modification anatasc derived from trans- 
mission spectroscopy. The values for A > 10 /xm are from Posch et al. (2003). 
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(Bradley 2003), and on planetary surfaces (e.g. Hartmann 2005). 

Olivine has the general sum formula (Mg,Fe)2Si04 with its Mg end member forsterite 
and the rarer Fe end member fayalite. On Earth, olivine is the most abundant material 
in the upper lithosphere. In general, more Mg-rich and less Fe-rich olivines are found. 
It crystallizes in orthorhombic symmetry with isolated [Si04]'*^ tcdrahedra surrounded 
by metal cations (Fe^+, Mg^+), each in coordination to six oxygen atoms. There are 
two different types of these coordinations, both slightly distorted from the octahedral 
symmetry. One is centrosymmetric and elongated along one of the 0-0 axes (cation 
sites in this coordination are called 'Ml'), and the other is non-centrosymmetric and 
irregular ('M2'; cf. Burns 1993). 

Olivine crystals usually have a greenish to yellowish color, which can turn into brown 
when it contains traces of Fe^+ ions. The greenish coloring is caused by a crystal field 
band of Fe around 1 /xm wavelength, which leads to an increased opacity for red light 
and an increased transmittance of yellow and green light. The brownish coloring comes 
from charge transfer processes toward the Fe^"*" ions leading to absorption increasing 
throughout the visible wavelength range towards the ultraviolet. This kind of absorption 
has been detected in the central parts of our Sri Lanka olivine sample (see below), where 
cracks may have given rise to oxidation of Fe^"*" in interface regions. 

In the present paper, we compare near-UV to MIR absorption spectra of fc(A) of San 
Carlos olivine single crystals with spectra of an Fe-rieh variety of olivine obtained from 
Ratnapura (Sri Lanka) (Fe contents: 8% and 20%, respectively see Table 3). Since 
olivine crystals have non-eubie symmetry, the polarization of the incoming radiation 
relative to the three crystallographic axes x, y, and z plays an important role in the 
resulting spectra. Several platelets of the San Carlos olivine have been prepared to have 
a surface perpendicular to one of these axes. Two of them, denoted as C and Dl, have 
been chosen to represent the properties of oriented olivine in polarized light due to the 
quality of their spectra. Sample C has been cut along the y-z plane, sample Dl along 
the x-y plane, so all three crystallographic axes were available for the measurements. 
The thicknesses are 1063 /itm for sample C and 1115 fxm for sample Dl. For the Sri 
Lanka olivine, only measurements of light polarized along the y- and ^-axes {E\\y and 
E\\z) were possible, since only one platelet (thickness: 300 iim) could be prepared. The 
sample is shown in Fig. 8. For these measurements we chose an area of about 1mm 
diameter without cracks, which can be seen in the figure. This area is also free of 
brownish colouring which occurs along the major cracks in the middle of the sample. 
The surface of this platelet is also not perfectly oriented perpendicular to the x-axis, 
because an inclination of approximately 20° could be measured by X-ray diffraction. 

The most important differences in the absorption spectra between the different ori- 
ented samples of San Carlos olivine are found in the bands eauscid by transitions of the 
metal cations, located around 0.7 1.5 /im. While E\\y and E\\z polarized light creates 
broad structures with peaks around 1.056 /xm and 1.108 yum, light polarized along the 
X-axis creates a more prominent band around 1.073 /im. Shoulders appear in all three 
polarizations on both sides of the main peak around 0.87 /xm and between 1.25 and 1.35 
/im. The value of k in the maxima of this band reaches ^1.8x10-"^ for both E\\y and 
E\\z, while k reaches a value of ^4x10^'' for E\\x. The corresponding values for the Sri 
Lanka olivine are given in Table 4. The k values of this band are increased compared to 
the measurements of the San Carlos olivine in the same polarization (in its maxima by a 
factor of «2) and also the particular positions of band and shoulders are shifted to longer 
wavelengths, except for the short- wavelength shoulder that shifts to shorter wavelengths. 

The continuum k values are in general relatively independent of the polarization and 
increases constantly towards longer wavelengths (see Fig. 9). Toward shorter wavelengths 
we cannot see an absorption increase, which indicates noticeable charge transfer pro- 
cesses. The sharp increase in k at around 0.4 /xm, which becomes more pronounced with 
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Figure 8: Platelet of the Fe-rich Sri Lanka sample used for transmission measurements 
under polarized light. The sample has a width of 4 mm and a length of 6 mm. The black 
circle features the spot of the measurements. 



Table 4: Specific positions and k values for the main peaks ('mp') and shoulders ('sh F 
and 'sh IF) of the crystal field band caused by Fe ions of the San Carlos (SC) and Sri 
Lanka olivines (SL), respectively. 
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^max 
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SC 
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1.074 
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1.232 


3.02 


E\\y 


0.882 


1.11 


1.057 


1.78 


1.234 


1.43 


E\\z 


0.858 


0.83 


1.114 


1.65 


1.323 


1.12 


SL 














E\\y 


0.838 


2.08 


1.069 


3.77 


1.255 


3.04 


E\\z 


0.831 


1.63 


1.14 


3.33 


1.378 


1.76 



increasing Fe content, is due to the onset of intrinsic transitions of the Fe^"*" ions. The 
continuum k value is independent of the Fe content of the olivine, which is remarkable. 

Towards longer and towards shorter wavelengths from the crystal field band at around 
1 /xm, several very weak bands appear. The bands at shorter wavelengths are compara- 
tively sharp and appear independently of the polarization at 0.65, 0.49, 0.47, and 0.45 
/im. These bands come from the spin forbidden transitions in the Fe^+ ions (Burns 1993). 
In contrast to that, the structures at longer wavelengths show a polarization dependence. 
For the Sri Lanka sample in E\\y polarization, a broad and relatively strong band can 
be seen at 3 /im, which is due to the stretching vibrations of OH groups. For the E\\z 
polarization this band is much sharper, similar to the corresponding band in the rutile. 
Around the 5.5 /im wavelength, another band is located, which is likely an overtone of 
the 10 /im stretching vibration band of the Si04-tetrahedron. 
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Figure 9: Spectra of oriented San Carlos olivine in polarized light. The spectra were 
measured from samples C and Dl (see text for further details). 
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Figure 10: Spectra of oriented San Carlos and Sri Lanka olivines in polarized light for 
E\\y. For the San Carlos olivine spectra of sample Dl are shown. 
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Figure 11: Spectra of oriented San Carlos and Sri Lanka olivines in polarized light for 
E\\z. For the San Carlos olivine, spectra of sample C are shown. 

4 The influence of k on the temperature of small grains 



In the present section, we apply our results for the magnitude of fc(A) to the calculation 
of the radiative equilibrium temperature of dust grains, which is an important parameter 
characterizing the conditions in dust-forming regions like circumstellar shells, and it 
indeed depends on absorption and emission properties. In the case of a thermodynamical 
equilibrium between the stellar radiation field (assuming a star with radius i?* with an 
effective temperature T*) and the dust grains in an optically thin circumstellar shell, it 
is possible to derive the dust temperature for a given distance from the star from the 
energy balance. The radiative energy absorbed by an individual dust grain is then given 
by 

-E'abs = -Drad h , (7) 

where 

h= 7ra2gabs(A)7rB(A,T*)dA, (8) 

^0 

and Dradi the radial dilution factor of the stellar radiation field, is defined as 



i'rad = 2[1 - 




Here, i?* denotes the effective stellar radius, and i?d is the distance of a dust grain 
from the star's center. Equation (9) is based on the so-called Lucy-approximation (Lucy 
1971), according to which the star is not considered as a point source, but its spatial 
extension is taken into account, which is important for small distances of the dust grains 
from the star. Already for a distance of 3 stellar radii, the dilution factor -Drad may be 
simplified to Rl/Rl- 

The radiative energy emitted by a single (spherical) dust grain with the radius a is 
given by 
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Eem = I2= 47ra2g^bs(A)7rB(A,Td)rfA. 
Jo 



(10) 



Numerically, the integrals li and I2 cannot (and don't need to) be calculated from zero 
to infinity as equations (8) and (10) suggest. It is sufficient to integrate over those wave- 
lengths where the star emits a signifieant fraction of it total radiative energy. On the 
other hand, the energy balance of the dust cannot be calculated correctly if optical con- 
stants (and hence Qabs- values) are lacking even at NIR wavelegths where the maximum 
of the stellar radiation field is located. 

The combination of eqs. (8) (10) leads to eq. (11), which allows calculating the de- 
pendence of the dust temperature on the distance from the star R^: 

i?d(Td) = . ^* , • (11) 

^1 - ti - -tt]' 

For the simplified dilution factor E^/R\, the previous equation can also be simplified, 
namely to 




Rd{Td) = \ -r R*- (12) 



For very small dust grains (with sizes amounting to 0.01 /im or less), the relation be- 
tween Td and Rd becomes independent of the grain size (which enters into the absorption 
efficiency Qa6s(A), in both the integrals li and I2; see Kriigel 2003). In this case, the 
infiuence of the optical constants on QabsW becomes most decisive for the dust temper- 
ature. Again in the small particle limit (also called Rayleigh limit), the relation between 
Qabs and the optical constants n and k is given by 

, , 27ra Qnk . 

yabs(A)=4— ^^^_^2_^2)2+4„2fc2' (13) 

from which it follows that in those regions where n(A) is approximately constant, such as 
in the NIR region for the oxide species discussed in this paper, A;(A) becomes the decisive 
quantity for Qabs (A), hence decisive for the dust temperature as well. 

More precisely, it is the magnitude of fc(A) and Qabs(A) in the visual and NIR relative 
to the respective MIR/FIR-values that determines the dust temperature. Consequently, 
dust species that have small kiX), corresponding to their high transparency in the optical 
and NIR - such as transition metal free oxides and silicates - will be much less heated 
by the stellar radiation than absorbing grain species. This can be seen most clearly from 
Fig. 12. For MgO, the lack of absorption mechanisms operating at short wavelegths 
leads to inefficient radiative heating, such that the dust temperature T^ drops below 
500 K already at two stellar radii. For Mgo.5Feo.5O, on the other hand, T^ amounts to 
2000 K close to 2R*. Furthermore, Td(R(j) drops much more steeply for MgO than for 
Mgo.5Feo.5O. 

The difii'erence between MgO and Mgo.5Feo.5O with respect to NIR absorbance and 
radiative equilibrium temperature is obviously an extreme case. Smaller, but still sig- 
nificant differences in T(i(r) show up for the other dust species discussed in the present 
paper, e.g. for spinel (see Fig. 13). 

As a consequence of the differences in the absorption efficiency factors calculated 
according to eq. (13) for small spherical grains - see insert in Fig. 13 - a difference in 
Td results. In this case, it amounts to about 25% at 2R, (the precise values are 850 K 
vs. 625 K at 2R*). While Qlabs calculated according to eq. (13) for small spherical grains 
(Fig. 13 insert) increases linearly with k, the increase in T^ in this case appears to be 
relatively weak, the gradient of T^ with R remains steep for the Cr-containing spinel. 
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Figure 12: Comparison of the radiative equilibrium temperature of MgO (characterized 
by extremely low fc(A) values in the NIR region) with the equilibrium temperature of 
Mgo.5Feo.5O (its fc(A) value is more than seven orders of magnitude larger in the NIR). 
The stellar effective temperature has been set to 3000 K. 



Notwithstanding, it is important to take the effect of kyis/NiR into account, since 
when it comes to condensation or evaporation, a few degrees Kelvin may decide about 
the existence or nonexistence of grains in a particular zone around a dust-forming star. 
Furthermore, there is still another parameter that linearly increases with k in the small 
particle limit, namely the radiation pressure efficiency factor Qpr- For particles that are 
smaller than the wavelength, Qpr « Qa6s> and thus an increase in k leads, via eq. (13), 
to a corresponding increase in more efficiently a dust grain absorbs radiation 

around 1 /im, the larger the radiation pressure it experiences and the velocity it can reach 
in a radiation-pressure driven stellar wind. For a detailed discussion of the influence of 
absorption and scattering on the radiation pressure acting on dust grains, we refer to 
Hofner (2008). 



5 Conclusions 

We examined the near-infrared and visual absorption properties of several Stardust 
analogs. Our focus has been on materials for which extremely low k values - or no 
k values at all - can be found in the literature for the above-mentioned wavelength 
range: spinel, rutile, anatase, and olivine. 

The examined samples were either natural terrestrial ones or synthetic crystals with 
a defined amount of absorption-enhancing impurities. In all cases, the impurity contents 
were such as to be compatible with cosmic elemental abundances, and in all cases, the 
derived k values were 1-2 orders of magnitude higher than the limit of the 'transparency 
regime' defined at the beginning of this paper (Eq. (1)). 

We suggest that even for cosmic dust grains which have smaller impurity contents 
than our analogs, k values below the limit of transparency will rarely be reached. In 
cases where no appropriate k data are available, the limit of the 'transparency regime' 
(fc w 10~^ at 1 yum and k w 10~^ at 10 /um) may serve as a heuristic lower limit for the 
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Figure 13: The radiative equilibrium temperature of natural (terrestrial) spinel dust con- 
taining impurities compared to the temperature of impurity-free Mg-Al-spinol according 
to Palik (1985-98). The stellar effective temperature has been set to 4000 K. 



absorption of natural materials including Stardust grains, although certain circumstances 
can produce very transparent natural minerals. 

Even though the 'real' k values of cosmic dust grains still has to be determinated by 
further investigations (e.g. on presolar grains), the above constraint seems more realistic 
than adopting k values dropping to 10~^ and below, which can be found in the solid-state 
physics literature and which mostly refer to synthetic, impurity-frcc materials. 

As an example, we showed that using k data for spinel taken from Palik (1985-98) can 
lead to underestimating the equilibrium temperature reached by this oxygen-rich dust 
species in a circumstellar radiation field. More realistic results can be obtained if optical 
constants with k values above the transparency limit are used. The same holds true for 
olivine, which is a much more abundant component of circumstellar dust. 

The magnitude of absorption also has an important infiuence on other quantities 
characterizing an expanding circumstellar shell, such as the radiation pressure. This is, 
however, beyond the scope of the present paper. 
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